Can we trust the temperature - misfit strain thin film phase diagrams? 
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It was demonstrated that phenomenological approach based on Landau thermodynamic theory 
widely used for modeling of ferroelectrics can lead to a substantial inaccuracy of description of ther- 
modynamics of ferroelectric thin films if non-linear electrostriction and non-linear elastic compliance 
are not taken into account. In this paper we demonstrate this with an example of the tempera- 
ture - misfit strain phase diagrams of a single-domain (OOl)-oriented BaTi03 thin film. Non-linear 
electrostriction and non-linear compliance tensors were calculated using first principles methods. 
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Strain engineering is a modern strategy to control and enhance material properties. It represents a technique which 
deals with semiconductor and ferroelectric thin films strained on a substrate. Such strain can be tensile or compressive 
and occurs due to a difference between lattice parameters of the film and that of the underlying substrate. In strain 
engineering, a strain appeared in the film shifts the paraelectric-to-ferroelectric transition temperature and thus can 
change the properties of the material such as dielectric and piezoelectric constants, remanent polarization, or even 
can induce the ferroelectricity in a non- ferroelectric material [![. To describe the ferroelectric state of a strained 
film, it was found to be very convenient to use temperature - misfit strain phase diagrams. First introduced by 
Pertsev et al. Q such diagrams present polarization state depending on temperature and strain appeared in the 
film after deposition it on the substrate. They have been developed for many classical ferroelectric materials such as 
BaTiO 3 0-li, PbTi0 3 [3,l3, and SrTiO 3 0, etc. 

Theoretical study of temperature - misfit strain phase diagrams can provide strain engineering and relevant exper- 
iments with guide lines. For example, the first phase diagrams developed by Pertsev et al. 0, [S] shortly stimulated 
experiments which confirmed the theoretical predictions, demonstrating huge misfit strain effects on ferroelectric 
properties of BaTi03 ||. 

Calculations of such diagrams involve Landau theory or ab initio microscopic approach. The first principles cal- 
culations 6| are more and more potent nowadays but require a lot of computations and there still exists a variety 
of problems where phenomenological approach can be valuable. Phenomenological approach is based on Landau- 
Ginzburg-Devonshire thermodynamic theory which is widely used for phase field modeling for calculations of 
single-domain Q and poly-domain Q thin film diagrams. However, there is a considerable difference between the 
results obtained by first-principle calculations and by calculations based on Landau phenomenological theory 0, [f| . 
What is the reason of such disparity? There is no doubt that to have reliable description of a strained system based on 
thermodynamic theory one has to know precisely the coefficients of thermodynamic energy expansion which are not 
always easy to get experimentally. However, if we suppose that we know these coefficients enough precise, can then 
the thermodynamic description of a strained system be trusted? This question was already raised |9| in the context 
of PbTi0 3 thin film of (lll)-orientation and it was suggested that a classical set of thermodynamic coefficients not 
including non-linear electrostriction can be insufficient for reliable description of a strained ferroelectric film. The goal 
of this paper is to clarify this issue to show that non-linear electrostriction and non-linear mechanical compliance can 
substantially affect the thermodynamics of thin ferroelectric films. 

The problem with the standard description of strained ferroelectrics can be illustrated qualitatively by a simple 
scalar model. Let us describe a ferroelectic with a Gibbs thermodynamic potential energy expansion keeping only one 
component of polarization P and stress a: 



G=-P 2 + ^P 4 + IP 6 - -a 2 - QP 2 o - — P 2 a 2 - RP 4 a - -a 3 , (1) 
2 4 6 2 ^ 2 3' W 
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where "ordinary" part represents Gibbs energy commonly used to describe ferroelectric systems and "high-order" 
terms (M, R - high-order electrostrictive coefficients and N - non-linear compliance) which are customarily neglected. 
This neglect can be readily justified for a bulk material but the situation for a clamped system is different. We can 
show this considering ferroelectricity in a clamped system, i.e. where the strain e equals the misfit strain eo- To 
obtain an effective potential G(T,P,€q) = G + ea, T - temperature, for the clamped system, we eliminate stress a 
using mechanical equation of state 

dG 

e = - — = so- + Na 2 + QP 2 + RP 4 + MP 2 a, (2) 
do 



and write the G(T, P, e ): 



Electro-mechanical interactions in such system lead to renormalizations of a, (3 and 7: 
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f MQ 



R NQ 2 



(5) 



Q 2 Q Q 3 

7* = 7 - 3M\ + 6R- + 2N\ . (6) 
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Atomic order-of-magnitude estimates show that the corrections due to high-order coefficients A and B are expected 
to be small as compared to corrections eo^ f° r a and 2 — for (3 (except the cases where the low-order corrections 
are unusually small) whereas C correction may be of the same order of magnitude as 7 ((6]). This situation is similar 
the strong renormalization of j3 term (charge of sign) in BaTi03 films caused by electrostriction 0] . These estimates 
suggest that the high-order coefficients, which are customarily neglected in majority of problems, should be included 
into thermodynamic energy expansion as far as 7 is important for the description of a problem. 

To illustrate this phenomenon quantitatively, we demonstrate the influence of high-order terms on the temperature 
- misfit strain phase diagram for a single-domain (OOl)-oriented BaTiC>3 thin film. To plot the diagram, we use 
Landau-Ginsburg-Devonshire thermodynamic theory with boundary conditions applied to a 2D clamped film on the 
substrate. We start with the Gibbs expansion with respect to polarization Pj and stress o~i (Hereafter we use the 
Voigt matrix notation and assume summation over repeated indices.): 

G = ai P 2 + a i: P 2 P 2 + a l3k P 2 P 2 P 2 - "-fa^j - —^-,0,0, - Q 13 P 2 a 3 - ^P 2 a,a k - R ljk P 2 P 2 a k , (7) 

where d{, aij, and aijk are dielectric stiffness and higher-order stiffness coefficients at constant stress, Sij and Nij k are 
linear and nonlinear elastic compliances, Qij ordinary electrostriction, and Mij k and Rij k are high-order electrostric- 
tion tensors. To obtain effective potential G(Pi, T, eo) = G + e\a\ + 6202 + e606, the minima of which correspond to the 
ground state of the film, we apply mixed boundary conditions ( J^- = e , = eo, §^ = and er 3 — 0, 04 = 0, og = 0) 
[3| and eliminate ai, 02 and eg. Hereafter Cartesian coordinate system with the X3 axis perpendicular to the film- 
substrate interface is considered, eo = a|l ao a ° - biaxial parent misfit strain, where ao is the lattice parameter of the 
ferroelectric in the cubic phase extrapolated to temperature T and an is the in-plane lattice parameter of the film. 
The function G obtained depends on Pi, T, and eo. Minimizing G = G(Pi,T, eo) with respect to polarization p for 
each (T, eo) point one can build the whole phase diagram. 

To see what impact high-order interactions exert on the strained BaTi03 thin film, we should know the values of 
high-order electrostrictive and non-linear compliance tensors themselves, which are not yet available experimentally. 
Consequently we turned towards ab initio methods, namely, we used Vienna Ab-initio Simulation Package (VASP) [10( 
performing zero Kelvin Density Functional Theory (DFT) full relaxation calculations. All calculations were performed 
within the generalized-gradient approximation as implemented in VASP using the projector augmented-wave method 



for the electron-ion interactions We have used a 6x6x6 Monkhorst-Pack grid for k-point sampling [12j, and a 

plane- wave energy cut-off of 600eV. For full relaxation calculations, the threshold of Hellman-Feynman force was less 
than lmeV/A. Here we would like to underline that mechanical and electro-mechanical properties, like mechanical 
compliance and electrostriction in ([7]), are expected to be weakly temperature dependent, what justifies the use of 
zero Kelvin DFT results in finite temperature calculations. 

As the first step, using VASP we obtained all components of the 6th-rank high-order electrostrictive tensors mi jk and 
rij k as well as of non-linear compliance n,ij k of Hclmholtz thermodynamic function F written in terms of polarization 
P and strain a: 



F = b t P 2 + h 3 P 2 P 2 + h ]k P 2 P 2 P 2 + 'f<,<, + >.<,<:.. - qaPfei - ^Pfe 3 e k - r l]k P 2 P 2 e k . (8) 



Technically, because of working with DFT, it is more convenient for us to use F expansion than G. When the c^-, 
Tiijk, Qij, n2ij k , and r-yfe coefficients are known, one can find coefficients Sy, Nij k , Qij, Mij k , and Rij k for expansion G. 
The relations between coefficients of G and P, in view of their complexity, are presented in Supplemental Materials 
at [URL inserted by publisher] . 

Let us briefly demonstrate the way we found stiffness and riij k and electrostrictive qij, raij k , and 77^ coefficients. 
Using VASP we find <r(e) and P(e) dependencies, then we put them to mechanical state equation 



<?h = c h j6j + nh^e-jCk - q ih P 2 - m^hPfe-j - r i:jh P 2 P 2 , 



(9) 
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and determine the coefficients by fitting. It is possible to find Cy and tensors separately from others (&y, m^fc, 
and Tijk) if one uses the mechanical state equation (|9]) at zero polarization and model deformation to paraelectric 
cubic phase with VASP (keeping mmm symmetry of the structure) . Then the other coefficients can be found from 
analysis of the full Eq. @. Further for simplicity we demonstrate how to find qu, mm and rm tensors components. 
One uses Eq. §§§ for the case where only polarization Pi and strain t\ are nonzero: 

o-i = cuei ± nmei - quPf - mm^ti - r nl P?. (10) 

Using VASP full relaxation calculations and keeping 4mm tetragonal symmetry of the structure, we find Pi(ei) and 
ci(ei) dependences (polarization was calculated by atom displacements and Born charges which were calculated for 
the cubic phase) and substitute them into Eq. (fTO]) . Once we have polynomial functions of t\ on the both sides of Eq. 
(fl~0|) we can obtain qu, mm, and rm. By repetition of such a procedure for different symmetries of the structure and 
for different components of <7j and 6j it is possible to calculate all the components of the high-order electrostrictive 
tensors. 

The coefficients of ordinary and high-order electrostriction as well as linear and non-linear elastic compliances 
obtained with the first principles calculations are given in the Table fl] Here, we should notice that current DFT 



Mnr 12 ^] 




M, [lO" 12 ^] 




n,[iq- 23 ^} 


Sll 


4.25 ±0.01 


Qu 


0.11 ± 0.01 


Mm 


-5.6 ±0.6 


Rm 


0.5 ±0.1 




8± 1 


S12 


-1.14 ±0.01 


Q12 


0.04 ±0.01 


M112 


0.8 ±0.5 


R121 


< 0.2 


Nll2 


< 2 


S44 


8.45 ±0.01 


Qu 


0.04 ±0.01 


Mx22 


1.6 ±0.6 


Rll2 


-0.8 ±0.1 


N123 


< 2 




M 123 


< 0.5 


R123 


-0.3 ±0.2 


iVl44 


< 2 


Q from [3],[f£] 




M144 


< 1.6 


-R144 


< 0.2 


NlS5 


< 2 


Qn 


0.11 


Qu 


0.15 


M155 


4.4 ± 1.4 


#156 


< 0.2 


N456 


< 4 


Ql2 


-0.043 


Q12 


-0.03 


M414 


< 1 










Qu 


0.059 


Q44 


0.04 


M 42 4 


< 1 


















M 45 6 


< 1 











TABLE I. Some material parameters of BaTiC>3 obtained from ab initio calculations, s^, Nijk, Qij, My*,, and Rijk are elastic 
compliance, nonlinear elastic compliance, electrostrictive tensor, and high-order electrostrictive tensors, respectively. 

methods give relatively moderate precision for high-order coefficients as it is clear from the error bars in Table HI 

As it was said above, discussing the scalar model, these high-order coefficients lead to changes of tensors of coeffi- 
cients dijk — S- a*j k of P 6 -terms while passing from a bulk crystal to the thin film (G —> G). Table [TTI demonstrate the 
size of this effect for BaTiOa (for detailed information on a*j k renormalization see Supplemental Materials at [URL 
inserted by publisher]) where the components of the original and renormalized tensors as well as the renormolizing 
corrections are given. "Corrections" represent an addition to cnjk, for example 0333 = am ± Aag 33 . 



Bulk 


Correction 


Film 


am = 8 


Aa* n = 4 
Aa^ 33 = -20 


a m = 12 ±3 
a^ 33 = -12 ±3 


au2 = 5 


Ao*i2 = -31 
Aa* 113 = -27 
Aa* 33 = -51 


alia = -26 ± 4 
ali3 = -22 ± 5 
ai33 = -46 ± 5 


ai23 = 5 


Aim = -62 


ai23 = -57 ± 20 



TABLE II. (Column 1) coefficients a ijk of P 6 -terms in 10 9 ^p at 300K of expansion G for bulk and mechanically free BaTiOs; 
(Column 2) corrections for corresponding a*j k ; (Column 3) a*j k of G of a (001) clamped film of BaTiOs. 

Inspection of this table shows that the renormarization is strong, confirming the conclusion drawn above from the 
order-of- magnitude estimates. Another striking feature of the data from the table is negative sign of some of a*- k , 
making the renormalized theory formally unstable. However, we found that physically the situation can be treated as 
stable. The point is that the potential G of thin film still has local minima due to the positive sign of renormalized 
coefficients for P 4 -terms of G (as in Ref. Q), therefore the system is locally stable. As for the global stability, it can 
be restored by adding P 8 -terms G expansion 13]. We do not incorporate these terms in our consideration, however 
we believe that this will not essentially affect the positions of the local minima in view of very high power of the 
terms. 

Figures QJa) and HJb) show how original phase diagram built with the coefficients from [3J changes when supple- 
mented with the high-order coefficients from our first principles calculations. Taking into account the high-order 
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FIG. 1. Temperature - misfit strain phase diagrams of a single-domain (OOl)-oriented BaTiC>3. eo - biaxial parent misfit strain, 
T - temperature. The phases are denoted as (i) the c-phase (yellow), polarization is out-of-plane (Pi = P2 = 0, P3 7^ 0); (ii) the 
aa-phase (red), where polarization is in-plane (Pi 7^ 0, P2 7^ 0, P3 =0); (iii) the ac-phase (brown) where Pi 7^ 0, P2 = 0, P3 7^ 0; 
(iv) r-phase (red), where all components of polarization are non-zero; and (v) paraelectric phrase (blue), polarization is 0. 
(a) Original Pertsev's diagrams built with coefficients from 3:]. (b) Developed with coefficients set from [j| appended with 
high-order Mijkimn, Rijkimn and Nijkimn coefficients. For (c) and (d) elastic and electromechanical couplings were obtained 
from first principles calculations, (c) Calculated with "traditional" electrostrictive coefficients used for description of all 
electromechanical couplings, no high-order terms, (d) Calculated with the full set of electromechanical coefficients. The 
hatched regions demonstrate the shift of the transition lines within the error bars of the coefficients. 



coefficients strongly changes the diagram: phase transition lines become shifted and ac-phase is disappeared like in 
6]. In the case of BaTiOa Mijk and AT^j. do not give large contributions while Rijk plays the major rule. 

An interesting fact we would like to mention is the positive sign of Q12 electrostrictive coefficient we obtained. At 
the same time, if we follow the "traditional" way of calculating electrostriction coefficients as Qi2 ad = ttt, where Ps 
is a value of spontaneous polarization and es is a value of spontaneous strain in tetragonal phase, our "traditional 
value" of <5i r 2 ad is negative and close to those previously reported 0,0] (see Table HJ. The implication is that the high- 
order electrostrictive coefficient -R112 was found anomalously large so that, when polarization and strain reach their 
spontaneous values, i?n2 strongly influences the "traditional" electrostriction which actually includes contributions 
of Q12 and i?n2: 

n trad £ S1 Ql2-P52 + Rll2Psj n D p 2 n \ 

^i2 = TT2 = 772 = via + ^ii2Ps 2 - I 11 ) 

PS2 PS2 

We double checked the positive sign of Q12, performing frozen phonon calculations to obtain a positive Q12 about 
0.03^. 

In view that electrostrictive coefficients were found different from those used in Q it's instructive to plot the 
phase diagram where all elastic and electro-mechanical coefficients are taken from ah initio calculations. Figures 
0Jc) and Old) show two diagrams based on ab initio calculations. Figure [TJc) is calculated within the standard 
framework (where all electromechanical effects are described with conventional electrostriction coefficients) using the 
" traditional" scheme for the treatment of our VASP data. Figure [ljd) is developed using the full set of coefficients 
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including high-order terms obtained from the same data. The hatched regions demonstrate the shift of the transition 
lines within the error bars of the coefficients, in the hatched region between c-phase and r-phase a new ac-phase can 
appear similar to Fig. [TJa). It is seen that the diagram is extremely sensitive to the variation of the high-order 
coefficients which have relatively low precision as it can be seen from Table HI 

Our analysis shows that an adequate Landau theory treatment of thermodynamics of ferroelectric thin films (single- 
domain and poly-domain thin films calculations, and phase field modeling) requires taking into account high-order 
electromechanical couplings and non-linear elasticity. In view of this finding, we believe that an experimental evalu- 
ation of high-order electromechanical couplings in ferroelectrics seems to be a tack of primarily importance. 
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SUPPLEMENTARY MATERIALS 



Knowing the coefficients of Helmholtz thermodynamic function F written in terms of polarization Pj and strain 

F = h.l'f + h.J'rl'j + b ijk PfPjPl + '-f,,,, + ^eiCjCk - QijPfe - ^fl'f'ji, - r, jL !'?!>;,,,. (12) 

it is possible to find the corresponding coefficients of Gibbs energy expansion G with respect to polarization P, and 
stress <ji 



G = a t Pf + ayPfP? + a ljk P?PfPt - -fa^ - -fL<j l( T 3 <j k - Q, ,,/>>., - —^I'f^ - R ijk PfPfa k . (13) 



~2 

One performs Legendre transformation 



G = P-^e l( r l (14) 



i=l 

where ti(Ji represents a work required to maintain a constant strain. Then using mechanical equation of state 

dF 

(Xi = -—,i = 1..6 (15) 



one eliminates stresses 6j between (fT4| and (fl5j) : solutions of Eqs. ([15]) were expanded in series and only low order 
terms of the expansion were kept, then Ci were substituted into (I14[) . 

The transformations of the corresponding high-order electrostrictive coefficients are listed below ([T6l - 130|) . 

M (di +ci 2 ) 2 min 4ci2 (en + ci 2 ) mug ,, fi ^ 

111 " (cii-c 12 ) 2 ( Cll +2c 12 ) 2 (cii-c 12 ) 2 ( Cll +2 Cl2 )2 + \ V 

2c\ 2 m 12 i 2c 2 2 toi 23 



(en - C12) 2 (en + 2c i2 ) 2 (en - C12) 2 (en + 2c i2 ) 2 ' 

M _ C12 (en +ci 2 )m m (c 2 ! + C11C12 + 2c 2 2 ) mn2 

112 " (cn - da) 2 (c u + 2c 12 ) 2 + (c u - d 2 ) 2 (c n + 2c 12 ) 2 1 <] 

CnCi 2 mi 22 CiiCi2m 123 



(cn - C12) 2 (cn + 2ci 2 ) 2 (cn - C12) 2 (di + 2ci 2 ) 2 

w = c 2 2 m m 2ciic 12 m 112 

122 (cii-c 12 ) 2 (c 11 +2 Cl2 ) 2 ( C ii-c 12 ) 2 ( Cll +2 Cl 2) 2 + 1 J 

(c 2 x + 2did2 + 2c 2 2 ) mi22 2ci 2 (di + c i2 ) m i23 



(di - C12) 2 (di + 2da) 2 (cn - da) 2 (di + 2ci 2 ) 2 : 



c\ 2 miii 2c 11 c 12 m 112 

(di - C12) 2 (di + 2d 2 ) 2 (cn - C12) 2 (di + 2c i2 ) : 
2ci2 (di + C12) TO122 (c 2 i + 2did2 + 2c 2 2 ) m 123 



Ml23 " ~ : 727^7^2 ^ TCT^TX^YT - ( 19 ) 



(di - C12) 2 (di + 2ci 2 ) 


2 ' (cii-ci2) 2 ( Cl i + 2c 12 ) 2 


M 144 = 


mi44 
9 ' 




c 44 ^ 


M 155 = 


TO155 

9 ' 






(dl + d2) 


2Ci 2 m424 



(20) 
(21) 



1 (cn - C12) (cn 4- 2cia) C44 (cn 2 + C11C12 - 2ci2 2 ) c 4 4 ' 
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M _ Ci 2 m 4 l4 gi 11)1424 

424 (cn 2 + C11C12 - 2ci2 2 )c 44 (en 2 + C11C12 - 2ci2 2 )c 44 ' 

M 456 = (24) 
c 44 2 

R _ (cn + C12) mm (gngii + c i2 (gn - 2g i2 )) 2ci 2 toi 22 (ci 2 gn - cngigj , . 

111 ~ (cii- C i 2 ) 2 (cii+2c 12 ) 2 ( Cll - Cl2 ) 2 ( Cll +2c 12 ) 2 + { ' 

2ci 2 mi 23 (ciggii - gngig) 2wn 2 (c 2 ^ + c u ci2 (~2gn + gig) + 2c 2 2 (-gn + gig)) 
(cn - C12) 2 (cn + 2ci 2 ) 2 (cn -ci 2 ) 2 (cn +2c i2 ) 2 

(cu + C12) r m 2ci 2 rii 2 



(cn - C12) (cn + 2c i2 ) c n + C11C12 - 2cf 2 

R = "^i2 2 (cfign + 2c 2 2 gn + C11C12 (gn - 4gi 2 )) cnmm (-giggii +cngi 2 ) . . 
121 ~ 2( Cl i-ci 2 ) 2 (cii + 2ci 2 ) 2 2(c 1 i-ci 2 ) 2 (ci 1 +2 C i 2 ) 2 + 1 ° j 

™i 2 3 (-3ciiCi 2 gii - 2c 2 2 (g n - 2g 12 ) + cfigi 2 ) mn 2 (2c 2 2 gn + c n (g n + 2g 12 ) - c n c 12 (g n + 4g 12 )) 



2 



2 (cn - C12) 2 (cn + 2c 12 ) 2 2(cn -ci 2 ) 2 (cii + 2c 12 ) 

ci 2 m 4 i 4 g 44 cnm 424 g 44 cnri2i ci 2 ri 23 



2 (c^ + C11C12 - 2c 2 2 ) c 44 2c 2 1 c 44 + 2cnci 2 c 44 - 4c 2 2 c 44 cf 1 + C11C12 - 2c 2 2 c 2 x + CuCi 2 - 2c 



12 



R _ wn 2 (c n gn + 2c 2 2 gn + cnci 2 (gn - 4g i2 )) _ c 12 m m (gngii + gig (gn - 2g 12 )) 
112 " (cii-ci 2 ) 2 ( C ii+2 C i 2 ) 2 (cii-c 12 ) 2 ( C i 1+ 2 C i 2 ) 2 

cnmm (-Ci2gn + cugi 2 ) C11TO123 (-ci2gn + cngi 2 ) ci 2 r m Cnrm 



(cn - C12) 2 (cn + 2c i2 ) 2 (cn - C12) 2 (cn + 2c i2 ) 2 cf 1 + C11C12 - 2c 2 2 c\ x + cuci 2 - 2c 2 2 ' 

„ _ ^123 (ciign + cncig (2gn - 3gi 2 ) + 2c 2 2 (gn - g i2 )) ci 2 m m (ci 2 gn - cngig) _ 
123 " (cii- C i 2 ) 2 (cii + 2 C i 2 ) 2 (cn-c 12 ) 2 ( Cl i + 2 Cl2 ) 2 + 1 Sj 

ram (c 2 igi2 + C11C12 (-2gu + gig) + 2c 2 2 (-gn + gig)) mug (cfigig + 2c 2 2 gi2 - cncig (2gn + gig)) 

(cii-ci 2 ) 2 (cii + 2ci 2 ) 2 + (cn- C i 2 ) 2 ( C ii + 2c 12 ) 2 + 

(cn + Ci 2 )m 4 i 4 g 44 Ci 2 TO 424 g 44 2c 12 ri 2 i (c u + c 42 ) r 123 



2 (cn - C12) (cu + 2ci 2 ) C44 (c n +ciici 2 - 2c 2 2 )c 44 c\ x + C11C12 - 2cf 2 (cn - c i2 ) (cu + 2ci 2 ) 



TO414 (cngn + ci 2 (gn - 2gi 2 )) -ci 2 m 424 gn + CnTO 424 gi 2 m 444 g 44 m 456 g 44 r 444 
2 ( c ii + c n c i2 - 2c 2 2 ) c 44 (cn - C12) (cn + 2ci 2 ) C44 2c| 4 c 2 4 c 44 ' 



™4i4(cngi2 - giggii) m 424 (-2ci 2 gi 2 + cn (gn + gig)) mi 55 g 44 ri 55 

(Cll - C12) (Cn + 2C12) C44 (Cll - C12) (Cn + 2C12) C44 C44 c 44 
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As it was discussed in the main text of the paper, to obtain effective potential of a film, G, 

G = G + ei<7i + e 2 cT2 + (31) 

on applies mixed mechanical conditions 

( 9G_ _ , dG_ _ , dG_ _ n 
J 9<Ti ~~ fc °' da 2 — fc 0, das — O 

10-3 = 0,0-4 = 0,0-5 = 0. 



and eliminates cti, ct 2 and 05 between (|3Tj) and (1321) . High-order couplings described by tensors Mijk, Rijk, and Nijk 
renormalize a^-fc coefficients of P 6 -terms of G. The renormalization reads: 

„ M m (0 llSll -Q 12 s 12 ) 2 M 122 (Qf 2 sf 1 

anl " anl 2( Sll -si 2 ) 2 (*ii+si2) 2 2( Sll -.s 12 ) 2 ( S ii + S i 2 ) 2 1 J 

#112 (-Ql2Sll + Q11S11S12 + Ql2Sllgl2 ~ Qlisj 2 ) _ #111 (-Qllgjl + Ql2SllSl2 + Q11S11S12 - Q12S12) _ 

(sn - S12) 2 (sn + S12) 2 (su - S12) 2 (sn + S12) 2 

MU2 (-QilSllSl2 ~ Ql2SllSl2 + Q11Q12 (sfi + 5 2 2 )) 



(Sll ~ S12) 2 (Sn + S12) ; 

(Qn + Q12) 



-(Qii 2 (sii 2 A^iii + s n (iV m - 3A^ii2) S12 + N : 



111S12 



3(sn -si 2 ) 2 (sn +S12) 3 

Q12 2 (sn 2 iV m + s n (JViu - 3iV U2 ) s 12 + iVins 12 2 ) - 
QuQu (sn 2 (iV m - 3JV U3 ) + 4s u iV m si2 + (iVm - 3JV U2 ) S12 2 )), 



M 122 g 2 2 Mi2 3 Q 2 2 , 2Q 12 i?ii2 , 2Q 12 3 (iVm + 3iV 112 ) 

3 333 - a m ~ 7 ; 77 ~ 7 i vJ H 1 1 57 i 73 ' v- 44 ; 

(S11+S12K (sn + S12) sn+si2 3 (sn + si 2 ) d 



„ _ M122 (Qn s n ( s n - 2si 2 ) +Q 2 2 si2 (-2sn + S12) + 2QnQi2 (s? x - S11S12 + s\ 2 )) 

a ii2 — a ii2 -7 rjT ; 70 (35) 

2 (sn - si2) z (sn + S12) 2 

M112 (QllQl2 (s U - 4snSi2 + S 2 2 ) + Q 2 ! (s 2 x - S11S12 + s\ 2 ) + Q\ 2 {s\ x - S11S12 + s 2 2 )) 



(sn - S12) 2 (an + S12) 2 
Mi 55 Q| 4 (s 2 ! - s\ 2 ) 2 M 424 (Q11 + Q12) Q44 . Q44-R155 (s?i - s 



2^2 
12J 



2 (Sn - S12) 2 (Six + Sl2) 2 s\ A (Sll+Si2)s 44 (Sn - S12) 2 (Sn + S12) 2 S 44 

Mm (2QiiQi 2 siis| 4 (sii - s 12 ) + Ql 2 s ii s 44 ~ 2Qn s n s i2S 2 4 ~ 2QiiQi2SiiSi 2 s| 4 - 2(3 2 2 snsi 2 s| 4 + Q n s i2 s 44) 

2 (sn - S12) 2 (sn + S12) 2 s\ A 

-Rill (-Ql2S n g 2 4 + Qllg 2 lgl2sl 4 + Ql2SllS 2 2 sl 4 - Qnfxgflf) 

(sn - S12) 2 (sn + S12) 2 s\ A 
(i?i2i(-2QnSnS 44 - 2Qx2S?x s 44 + 2QiiSnSi2S 44 + 2(2i 2 s 2 1 si 2 s 44 + 2QnSns 2 2 s 44 + 2(2i 2 siis 2 2 s 44 

2Qiis 3 2 s 44 - 2<3i 2 Si 2 s 44 ))/ ((sn - 812) 2 (sn + S12) 2 s 44 ) 

#112 ( — QllS n s 44 + Ql2 s Xl s 12 s 44 + Qll s ll s 12 s 44 — Ql2Sl2 s 44) 



(Sll - S12) 2 (Sn + S12) 2 s\ 



-11 



{QU + (h2) (Q 4 4 2 (sn 2 - Si2 2 ) 2 A^155 + (gil 2 (sil 2 iVll2 " «ll(JViu + JVua)« 12 +7Vll2Sl2 2 ) + 



(sn - S12) 2 (sn + S12) 3 s 44 2 

Qi2 2 (sn 2 iVii2 - «ii(iVui + iVn 2 )si2 + iVii2Si2 2 )+QiiQi2(sii 2 (iV m + iVn 2 )- 4fluJViia*ia+(JViu + A r ii 2 )si 2 2 ))s 4 4 
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„ MuiQi2 (Qnsii - Q12S12) M123 (QnQi 2 Sn - Qnsnsi2 - Q?2 s n s i2 + QnQi2S? 2 ) ,„_,. 
a ii3 — a H2 / w ; 7 TT7 ; \~5 V-JOJ 

(Sn - S12) (Sn + Sl 2 ) ^ (Sn - S12) 2 (Su + S12) ^ 

M112 (guQl2£ll + Ql2 s ll ~ 2QuQi2£llSl2 - 2Qfgfll£l2 + QuQlgfjg + Ql 2 s\ 2 ) _ 

(sn - S12) 2 (su +S12) 2 
#123 (-2Qi2£n + 2Qns 2 1 si2 + 2Q 12 s 1Y s\ 2 - 2Qusf 2 ) _ 
(sn - S12) 2 (sn + si 2 ) 2 

#121 (-2Qnsf 1 + 2Qi 2 5 2 1 Sl2 + 2QiiSn5 2 2 - 2Q 12 sj 2 ) _ 
(sn - S12) 2 (sn + S12) 2 

#111 (~Ql2gll + Ql25n5l2 + Ql2gllgl2 ~ Qlgflg) #112 (-Ql2«ll + Ql2£ll£l2 + Ql2£ll£l2 ~ Qlgflg) _ 

(sn - S12) 2 (sn + S12) 2 (sn - S12) 2 (sn + S12) 2 

Mm (ggiigiafu (~3sn + s 12 ) + Q 2 x (s 2 ! + s 2 2 ) + Q\ 2 (3s 2 n - 2s n s 12 + sj 2 )) 

2 (sn - S12) 2 (an + S12) 2 

Qia 



r(4 QiiQi 2 (sn 2 iVii2 - s n (iV m + AT 112 )si2 + A r ii 2 si2 2 )+Qii 2 (sii 2 (-Win + -W112) 



(su - S12) 2 (sn + S12) 3 

4s u iV 112Sl2 +(i\r ul + iYii 2 )si2 2 )+Qi2 2 (sii 2 (A r iii + N U2 ~ 4:S n N n2 s 12 +(N 1 n + N 112 )s 12 2 )), 



* = M ni Q 2 2 MmQfg Mi 23 Qi 2 (Qii + Qi2) M122Q12 (2Qn + 3Qi 2 ) 

° 133 ° 112 2(,S1 1+Sl2 ) 2 (S11+.S12) 2 (S11+S12) 2 2(,si 1+Sl2 ) 2 

#112 (QllSll + Ql2Sll + QllSl2 + Ql2Sl2) . #121 (2Ql2Sll + 2Q12S12) #123 (2Ql2Sll + 2Q12S12) 



(S11 + S12) 2 (Sll + Si 2 ) 2 (S11+S12) 2 

gi2 2 (gii + Qi 2 )(iViii + 3iVi 1 2) 

(sii + S12) 3 



2Afii 2 Qi 2 (Qii +Q12) 2M111Q12 (-sn + 512) (gigfn ~ gufig) 
(S11+S12) 2 (su - S12) 2 (sn + S12) 2 

M123 (-4QllQl2*ll*12 + Qll («?l + S12) + Q? 2 + *?a)) 



(fill - S12) 2 (su + S12) 2 

2M i22 (Q11S11S12 + Q\ 2 (2an - si 2 ) si 2 - Q11Q12 (2s 2 x - S11S12 + s 2 2 )) MiuQl^ (s 2 4 - s\ 2 ) 



(sn - S12) 2 (su + S12) 2 2 (sn - S12) 2 (sii + S12) 2 sf 4 

2M424Q12Q44 2Q4 4 i?i44 (s 2 4 - S 2 2 ) 2 

(Sll + S12) S44 (sn - S12) 2 (Six + S12) 2 S44 
#123 (4QXXSXX (Six + SX2) (—Sll + £12) sl 4 + 4Qi2 (-Su - S12) SX2 ("Sn + Sx2 ) S44) 

(sxi - S12) 2 (sn + S12) 2 s| 4 

(#12l(8Qx2SxX (sn + S12) (-Su + S12) S 44 + 4Qu (-s u - S12) S i2 (-Su + S12) s 44 + 

4Qi2 (-su - S12) S12 (-su + S12) s 44 ))/((sn - S12) 2 (su + S12) 2 s 44 ) + 
2Q12-W111 (2QiiQi 2 sii 2 s 44 2 - 2Q 1 i 2 sii5i 2 s 44 2 - 2Qi 2 2 s 1 is 12 s 44 2 + 2QnQi 2 si2 2 s 44 2 ) 

(su - S12) 2 (su + S12) 3 s 44 2 

2Q12Q44 2 (sn 2 - S 12 2 ) 2 jVl55 | 
(SU - S12) 2 (Su + S12) 3 S44 2 
(2Qi 2 iVii2(2(3ii 2 Sii 2 S44 2 + 2QiiQi 2 Sii 2 S44 2 + 2Qi2 2 Sii 2 S 4 4 2 - 2Qii 2 SiiSi 2 S44 2 - 

8Q11Q12S11S12S44 2 - 2Q 12 2 siiSi2S44 2 + 2Q n 2 si2 2 S44 2 + 2QiiQi 2 si 2 2 s44 2 + 2Qi2 2 si2 2 s 4 4 2 ))/ 

((Su - S12) 2 (Su + S12) 3 S44 2 ) 

Note that in p2p solutions were expanded in series and only low order terms of the expansion were kept, then a\, a 2 
and (T 6 were substituted into (I31[) . therefore expressions (1331 - 138[) result from the expansion in series in (|32[) . 

Table Hill shows more detailed information for corrections Aa^ 11; Aa* 12 and AaJ 2 3 demonstrating partial contribu- 
tions from Mijk, Rijk and Nijk tensors. 
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8 Ao m , 


R, [fjr] 


8 Aan 2 , 


R, [frr] 


$ Aai 23 , 

[io 9 <^] 


a m k = 8 


«112 = 5 


a 123 k — 5 


Mm -5.6 ±0.6 


2.4 


Mm -5.6 ±0.6 


3.6 


Mm -5.6 ±0.6 


2.3 


M112 0.8 ±0.5 


-0.4 


Mna 0.8 ±0.5 


-1.3 


M112 0.8 ±0.5 


-1.0 


M122 1.6 ±0.6 


-0.2 


M122 1.6 ±0.6 


-1.5 


M122 1.6 ±0.6 


-2.7 


M123 < 0.5 j 





M123 < 0.5 





M123 < 0.5 


1.2 


M 144 < 1.6 





M 144 < 1.6 





M144 < 1.6 





M155 4.4 ±1.4 





M155 4.4 ±1.4 


-0.1 


Miss 4.4 ± 1.4 





M414 < 1 





M 4 14 < 1 





M 4 14 < 1 





M 4 24 < 1 





M 424 < 1 


-0.2 


M 424 < 1 


-0.1 


M 45 6 < 1 





M 456 < 1 





M 45 6 < 1 





#m 0.5 ±0.1 


14.7 


#m 0.5 ±0.1 


8.4 


#m 0.5 ±0.1 





R121 < 0.2 





R121 < 0.2 


-18.5 


#121 < 0.2 


-23.2 


#112 0.8 ±0.1 


-13.5 


Rii2 0.8 ±0.1 


-23.6 


#112 0.8 ±0.1 





#123 -0.3 ±0.2 





#123 -0.3 ±0.2 





#123 -0.3 ±0.2 


-35.3 


#144 < 0.2 





#144 < 0.2 





#144 < 0.2 


0.1 


#155 < 0.2 





#156 < 0.2 


< 0.5 


#155 < 0.2 





JV m 8 ± 1 


1.3 


Nm 8 ± 1 


0.3 


Nm 8 ± 1 


-0.3 


M12 < 2 


0.1 


N 112 < 2 


1.1 


N 112 < 2 


-1.1 


iVl23 < 2 





iVl23 < 2 





iVl23 < 2 





JVl44 < 2 





iVl44 < 2 





iVl44 < 2 





iV 155 < 2 





iVl55 < 2 


< 0.1 


iVl55 < 2 


< 0.1 


^456 < 4 





^456 < 4 





^456 < 4 





Total 


4.5 


Total 


-31.2 


Total 


-61.6 



TABLE III. Detailed information for corrections Aa m , Aa U2 and Aal 2 3- 8 means a contribution from a particular high-order 
electrostrictive coefficient. Values of a°^ k coefficients in 10 9 ^p for bulk BaTi0 3 at 300K are also give for comparison. 



